Energy functional for the three-level Lipkin model 
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We compute the energy functional of a three-level Lipkin model via a Legrendre transform and 
compare exact numerical results with analytical solutions obtained from the random phase approx- 
imation (RPA). Except for the region of the phase transition, the RPA solutions perform very well. 
We also study the case of three non-degenerate levels and again find that the RPA solution agrees 
well with the exact numerical result. For this case, the analytical results give us insight into the 
form of the energy functional in the presence of symmetry-breaking one-body potentials. 
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I. INTRODUCTION 

Nuclear density functional theory (DFT) is applicable 
for a description of ground-state properties throughout 
the nuclear chart 

[ESS, II- This theory is based on the 
theorem by Hohenberg and Kohn Q , and its practical ap- 
plications are based on the self-consistent mean-field the- 
ory with density-dependent energy functional [E EL II, U • 
Within nuclear DFT, the solution of the quantum nuclear 
many-body problem is relatively simple. The problem is, 
of course, with the construction of the energy functional. 

Most approaches to the nuclear energy-density func- 
tional are empirical in nature, and the inclusion of pairing 
properties [Tot fill. Il2l. Il3l. 113. Il5| is particularly important 
and challenging. There are only a few non-empirical ap- 
proaches to the nuclear energy-density functional. The 
density-matrix expansion pioneered by Negele and Vau- 
therin [ijj] expands the non-diagonal one-body density 
matrix along its diagonal, and the expectation value 
of a Hamiltonian thus leads to an energy-density func- 
tional. In another formal approach [17] . the energy- 
density functional is constructed as the Legendre trans- 
form of the ground-state energy as a functional of ex- 
ternal one-body potentials. This approach is closely re- 
lated to the Hohenberg-Kohn theorem [f|. However, its 
starting point, i.e. the ground-state energy of a quan- 
tum many-body system as a functional of any external 
potential, is only available for a few solvable or weakly 
interacting systems. Furnstahl and co-workers [H, Il9j 
followed this path and derived energy-density function- 
als for dilute Fermi gases with short-ranged interactions. 
For Fermi gases in the unitary regime, simple scaling ar- 
guments s ugg est the form of the energy density func- 
tional HIImSS,!!!. 

In nuclear physics, energy-density functional theory is 
a practical tool that is popular due to its computational 
simplicity and success. The universality of the functional, 
i.e. the possibility to study nucleons in external poten- 
tials, is seldom used. This distinguishes nuclear DFT 
from DFT in Coulomb systems and makes alternative 
formulations worth studying. For computational simplic- 
ity, we would like to maintain the framework of an energy 



functional. However, there is no need to focus on func- 
tional of the density. The description and interpretation 
of nuclear systems is often based on shell-model orbitals 
rather than densities. For instance, in nuclear structure, 
one is often more interested in the occupation of a given 
shell- model orbital or the isospin-dependence of the effec- 
tive single-particle energies than in the shape of the den- 
sity distribution. This language is natural for shell-model 
Hamiltonians that are based on single-particle orbitals. 
It is the purpose of this paper to study and construct the 
energy functional of such a system. 

Recently, the energy functional (in terms of occupa- 
tion numbers) was constructed for the pairing Hamilto- 
nian [24. In this paper we consider the three-level Lipkin 
model [25l. l2ol |27I . l28j . This is another solvable model 
that is relevant for nuclei. We will be able to gain insight 
into the formulations of energy functionals for two-body 
interactions that exhibit a continuous symmetry, and will 
study the effect of symmetry-breaking one-body poten- 
tials. 

This paper is organized as follows. In Sect. [H] we 
construct the energy functional of the three-level Lip- 
kin model with degenerate excited levels. We study the 
implications of symmetry-breaking one-body potentials 
in Sect. IIIH and conclude with our Summary. 



II. DEGENERATE ENERGY LEVELS 

We consider the three-level Lipkin model |2 5| . | 26l. [27], 
HH and follow the work by Hagino and Bertsch |2°|. This 



model consists of N fermions that are distributed over 
three levels (labeled as 0, 1, and 2) consisting of N de- 
generate, single-particle states each. The single-particle 
energies are £o = 0, and the two excited levels are degen- 
erate, E\ — £i = e. The Hamiltonian is 



H = 



e(hx 
V 



{K X K X + K 2 K 2 



K\K\ 



(1) 
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Here, 



N 

£4 

i=l 



a = 0,1,2 



is the number operator of level a, while 



N 



CQi 



a = 1,2 



(2) 



(3) 



transfer fermions from the level to the level a > 0. We 
assume e > for the spacing to the degenerate levels 1 
and 2. The operators c ai and c a i create and annihilate a 
fermion in state i of level a. The Hamiltonian is invari- 
ant under the simultaneous orthogonal transformation of 
the iV orbitals belonging to each level, and this symme- 
try facilitates the numerical solution of this problem [28| . 
Note that the Hamiltonian {1} is invariant under orthogo- 
nal transformations of states belonging to the degenerate 
levels 1 and 2, i.e. 



-u 



COS( 

- sin i 



sin i 

COS( 



(4) 



This symmetry will be broken by the Hartree-Fock (HF) 
solution (see below), and we will study the explicit break- 
ing of this symmetry by additional one-body terms of the 
Hamiltonian in Sect. Hill 

The energy functional is the Legendre transform 



Here 



F(n) = E(e) - en(s) . 



8E(e) 

n{£) = -BT 



(5) 



(6) 



is the ground-state occupation of levels 1 and 2, while 
E(e) is the ground-state energy as a function of the spac- 
ing e. The exact ground-state energy can be obtained 
numerically by diagonalizing a matrix of modest dimen- 
sions [28l | . and one can thus construct easily the exact 
energy functional. However, we are interested in gain- 
ing analytical insight into the problem. For this rea- 
son, we employ the analytical expression by Hagino and 
Bertsch [29[ that expresses the ground-state energy in 
terms of the HF energy and corrections due to the RPA. 
The HF result depends on the size of the parameter e/v 
where 



v = V(N-1) 



(7) 



is the effective strength of the two-body coupling V in 
the iV-body system. The spherical phase is found for 
e/v > 1, while one deals with a deformed phase for e/v < 
1. In the spherical (deformed) phase, the HF solution 
preserves (breaks) the symmetry ((4]) of the Hamiltonian. 
The results are 129] 



E 



for e > v 



(8) 
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FIG. 1: (color online) Comparison of the exact ground-state 
energy to the RPA approximation for a system of N = 20 
fermions. The RPA result shows a cusp at e = v = (N — 1) V , 
corresponding to the boundary between the spherical (e > v) 
and the deformed phase (e < v) 



E = — 



N + 4 
2N- 1 



N - I e 



for 



< v 



(9) 



Figure [T] compares the approximate ground-state ener- 
gies ((HJ and (J9j) with the exact result as a function of 
e/V for a system of N = 20 fermions. The approximate 
solution exhibits a cusp at the boundary e/v — 1 be- 
tween the spherical and deformed phases, but it agrees 
well with the exact result away from the phase boundary. 
Note that the cusp renders the RPA solution non-convex, 
and its Legendre transform is therefore not possible close 
to the phase boundary. Within the RPA formalism, the 
occupation number ^ attains non-physical negative val- 
ues close to the cusp. We will avoid this problem in the 
construction presented below. 

Let us start with the analytical construction of the 
energy functional ([5]) in the spherical phase. Here, the 
RPA solution is convex. We have 



dE 
~de~ 



- 1 



(10) 



This equation can easily be solved for e(n), and the func- 
tional thus reads 



F(n) 



\v\\fnlj 



(11) 



Thus, for small occupation numbers (i.e. the weak 
coupling limit v <C e), the functional is nonanalytical 
and exhibits a square root singularity. This is a rather 
generic feature of energy functionals and is also seen in 
the energy functional for the pairing Hamiltonian. 

We next turn to the deformed case (e < v) . The occu- 
pation number is 



n(e) 



2N- 1 



N — 1 e 
2 v 



VW 2 



(12) 
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FIG. 2: (color online) Comparison of the exact energy func- 
tional (solid line) and RPA-obtained functional (dashed line) 
for N — 20 fermions. Inset: The phase transition is clearly 
seen at the intersection of the deformed and spherical RPA 
solutions. 



FIG. 3: (color online) Ground-state energy as a function 
of occupation number obtained from the exact functional 
(solid line) and the RPA functional (dashed line) for N = 20 
fermions. The RPA solution fails at the phase transition. 



Equation (|T2f is difficult to invert analytically because of 
the term involving the square root. Note that this term 
also renders the occupation number negative as e ap- 
proaches v. We avoid this problem by expanding Eq. (fT2|) 
in powers of e/v, approximate e/v -C 1, and keep up to 
orders 0(e/v) in the resulting expression. The result is 



. . JV- 2n-l/2 

ein) = =-i> 

N-1 + V2 



(13) 



Insertion of this result and Eq. (J9j> into Eq. ([5]) yields the 
density functional in the deformed phase. 

Note that the energy functional is analytical in the 
strong coupling limit v 3> e. Again, this is a generic 
feature of energy functionals and was also exhibited in 
the pairing problem. 

Figure [2] compares the exact energy functional with the 
RPA solutions obtained in the spherical and deformed 
phase for N = 20. Note that the spherical (deformed) 
phase corresponds to sufficiently small (large) occupa- 
tion numbers. The inset shows that the RPA solution 
has a kink at the critical occupation number n f=a 0.7 
at the phase transition. Note also that the solution of 
the deformed phase is a good approximation over the en- 
tire range n > 0.7, although its derivation employed the 
approximation e -C v corresponding to n ss N/2, see 

Eq. mD. 

How severe is the presence of the kink at n sw 0.7 in the 
functional? To address this question, we employ the RPA 
functional in practical calculations. We add the one-body 
term en to the functional and numerically determine the 
occupation number 71 that minimizes the ground-state 
energy 



E(n) = F(n) 



en 



(14) 



The energy taken at the minimum is plotted as a func- 
tion of the minimizing occupation number n in Fig. [3] 
The result is also compared to the exact numerical result. 
We see that the energy functional, obtained via RPA, 



provides a good prediction of the ground state energy 
function of the three-level Lipkin model. At the phase 
transition, however, the functional predicts multiple val- 
ues for the ground state energy, due to the discontinuity 
in the derivative at the phase transition. Thus, for ex- 
ternal potentials e ~ v the energy functional does not 
provide good predictive power, and the relatively small 
kink visible in the inset of Fig. [5] is, of course, the reason 
for this shortcoming. 



III. NON-DEGENERATE ENERGY LEVELS 
AND SYMMETRY-BREAKING 

The two-body interaction and the one-body terms of 
the Hamiltonian ([1]) are invariant under the symmetry 
transformation ^ . In this section we study the breaking 
of this symmetry by a one-body potential, i.e. we lift the 
degeneracy of levels 1 and 2 in the Hamiltonian |T]) and 
modify its one-body term as follows 



e(hi + n 2 ) — > eihx + £2^2 • 



(15) 



The two-body interaction remains the same. Technically 
it is convenient to substitute e — > (e\ + £2)/% in Eq. ((T|) 
and to add the term 



£1 - £2 



(fix - n 2 ) 



(16) 



to the resulting expression. This approach is particularly 
convenient for the RPA calculation. 

The question of symmetry-breaking is an important 
one. The energy functional corresponds to the symmetry- 
preserving two-body interaction that is probed by a 
symmetry-breaking one-body potential. Thus, the func- 
tional depends on two variables (the occupation numbers 
of levels 1 and 2), and we are interested in the form of 
this functional and in its relation to the functional of 
the symmetry-preserving case. We thereby hope to get 
insights into energy functionals for atomic nuclei. The 
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nuclear shell model, for instance, employs a spherically 
symmetric two-body interaction, but deformed nuclei are 
often described within a deformed mean-field basis that 
breaks the rotational invariance through one-body terms. 

The numerical calculations are simple and require only 
a minor modification of the matrix elements correspond- 
ing to the one-body terms. For the analytical construc- 
tion of the functional, we have to perfom the HF and RPA 
calculations. We adapt the HF calculation of Ref. [29[ 
to the case of non-degenerate levels and obtain the HF- 
energy surface 



E(a,f3) = N sin 2 ot{e\ cos 2 (3 + e 2 sin 2 0) 



-vN sin 2 a cos 2 a 



(17) 



E(e 1je2 )/V 




20 



as a function of the two angles a and (3 of the HF trans- 
formation. In the case £\ = £2 of degenerate levels, the 
energy surface (fl7[) becomes independent of (3. 

In what follows, we assume £1 < £2. The minimum 
of the energy surface (jTTJ) occurs at a = (where j3 is 
arbitrary and can be chosen as zero) for e\ > v, and at 
cos 2a = E\jv and (3 = for E\ < v. Again, we refer to 
these cases as the weak-coupling regime and the strong- 
coupling regime, respectively. The Hartree-Fock energy 
thus becomes 



Ehf = 



, 



N 



IV, 
T £l 



J -e\lv 



£1 > v 
£1 < v 



(18) 



Thus, within the HF approximation, all particles stay 
in level in the weak-coupling regime, while the fcrmions 
occupy level and level 1 in the strong-coupling regime. 
The RPA equation takes the well-known form 



A B 

-B -A 



(19) 



The nonzero matrix elements of the 2x2 matrices A and 
B are 

3 

An = £1 cos 2a + —v sin 2 2a , 

• 2 i> . 2 „ 
^5-22 = £2 — £1 sin a + — sm 2a , 

Bu = — w(cos 4 a + sin 4 a) , 
B22 = —v cos 2 a . 

Solving for the eigcnfrcqucncies u yields 



£1 > v 



2(v 2 - e{) , £1 < v 

e\ — v 2 , £1 > v 

(£ 2 - £i)(£ 2 + V) , £1 < V 



(20) 
(21) 



FIG. 4: (color online) Ground-state energy E(ei,E2) as a 
function of the single-particle energies si t 2 from the RPA 
(mesh of dashed lines) and the exact result (dots) for TV = 10 
fermions in the three-level Lipkin model. The flat region 
corresponds to the weak-coupling regime, and the cusp of 
the RPA results marks the boundary to the strong-coupling 
regime. 



and the result thus becomes 
2 

E = 2^(\/ £2 J~ v2 ~ £ i) f ° r £1>V > (23) 
N + A 1 



I -v + -V£2 - EiVv + £2 



v 2 - e\ 2N+1 



2 4 
£2 _ N-l ej 
2 4 v 



£1 



for Ei < v . 



(24) 



In the weak-coupling regime, the energy (123j) is very sim- 
ply related to the energy © for the spherical phase, while 
the energy of the strong-coupling regime differs mainly 
from its deformed counterpart ((9|) through the term that 
is nonanalytical in the level splitting £2 — £1. It seems 
that the non-analyticity is an artifact of the RPA approx- 
imation. 

Fig. H] compares this analytical expressions (f2l?|) and 
(f2~4")l with the exact energy for a system of N = 10 
fermions. Again, the RPA fails at the boundary between 
the regime of strong and weak coupling, and the energy 
is not a convex function of the single-particle energies. 

The energy functional results from the Legendre trans- 
form 



F(n 1 ,n 2 ) = E(ei,e 2 ) 



2 



(25) 



The total energy is 



/* — E] 



HF 




TrA 



(22) 



and it is understood that ej — £j(ni, TI2) results from the 
inversion of the occupation numbers 



dE 



(26) 
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We start with the weak-coupling regime Si > v and 
find 



(27) 



Thus, the occupation number of the two energy levels 
depends only on the energy of the level itself, and the 
strength of the two-body interaction. Note that the oc- 
cupation numbers are symmetric under exchange of level 
1 and level 2. The inversion of Eq. (f2"T)) is straightfor- 
ward, and we obtain the energy functional 



2 

F(m,n 2 ) = - \v\^2yn 3 (n 3 + 1) 



(28) 



Note that the functional (|28|) is symmetric under ex- 
change of level 1 and level 2, and that it is simply the 
sum of two functionals. Its form could have almost been 
guessed from the functional (|TT|) for degenerate levels. 

We turn to the strong-coupling regime e\ < v which 
once again proves more difficult. The occupation num- 
bers 



ni 



1>2 



2N + 1 x TV -lei 



El 



v y/2(v 2 - e{) 



(x + x^ 1 -2) 



(29) 



(30) 



do not decouple as in the weak-coupling limit. Here, we 
employed the shorthand 



£2 



£2 - £l 

We solve Eq. (|30p for x and obtain 



x{n 2 ) = 2n 2 + l + 2y/n 2 {n 2 + 1) 



(31) 



(32) 



We insert this result into Eq. (|2"9"|) , approximate e% <C 
v and solve for e\. This expansion again renders the 
resulting functional convex, and we obtain 



. N -2ni + l/2- x/2 
£l(ni ' n2)= N-1 + V2 V 

Finally, we insert this result into Eq. (|3"Tj) and find 



e 2 (n-i,n 2 ) 



V + X E\ 

x 2 — 1 



(33) 



(34) 



The insertion of the expressions (|33|) and (|34[1 into the 
functional ([23)) thus yields the desired expression. Fig- 
ure shows the resulting functional and compares it to 
the exact solution. The agreement between the exact re- 
sult and the RPA result is quite satisfactory. Note that 
the boundary between the regimes of strong and weak 
coupling is a one-dimensional line in the region of small 
occupation numbers 774.2 < 1- 



F(n 1 ,n 2 



)/V 
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FIG. 5: (color online) Energy functional F{n\,ri2) of the 
three-level Lipkin model with non-degenerate single-particle 
levels for N = 10 fermions. The exact result is shown as dots, 
while the RPA result is shown as a mesh and the correspond- 
ing contour. 



IV. SUMMARY 



We constructed the energy functional for the three- 
level Lipkin model with degenerate and non-degenerate 
excited levels. Our analytical results are based on RPA 
calculations and subsequent Legendre transformations. 
They agree well with exact numerical results in the limits 
of strong and weak coupling, respectively, but the RPA 
fails at the boundary between both regimes. In partic- 
ular, the RPA energy is not a convex function of the 
single-particle energies in the strong-coupling regime. 

In the case of non-degenerate excited levels, a conti- 
nous symmetry is broken by single-particle terms of the 
Hamiltonian while the two-body interaction alone pre- 
serves this symmetry. In the weak-coupling regime, the 
energy functional obtained from non-degenerate levels is 
simply related to the functional obtained for degener- 
ate excited levels. Both exhibit a square-root singularity 
for small occupation numbers. The relationship is more 
complicated in the strong-coupling regime. However, the 
derived results should be useful in the construction of an 
occupation-number based energy functional for nuclear 
masses. Recall that the occupation number-based mass 
formula by Duflo and Zuker [301 * s superior to mass table 
calculations that employ nuclear energy-density function- 
als [H . This gives prospect to the development of a global 
nuclear energy functional based on shell model occupa- 
tions. 
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LLC (Oak Ridge National Laboratory). 
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